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Abstract. We have measured the damped motion of a trapped Bose-Einstein 
condensate, oscillating with respect to a thermal cloud. The cigar-shaped trapping 
potential provides enough transverse confinement that the dynamics of the system are 
intermediate between three-dimensional and one-dimensional. We find that oscillations 
persist for longer than expected for a three-dimensional gas. We attribute this to the 
suppressed occupation of transverse momentum states, which are essential for damping. 
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1. Introduction 

Trapped, ultracold gases offer a versatile way to investigate quantum many-body 
physics. Well-isolated from their surroundings, they can be controlled to cover a wide 
parameter space, giving access to regimes beyond the reach of other condensed matter 
experiments [TJ. Confinement reduces the dimensionality of a gas when the atoms have 
insufficient energy to reach excited quantum levels. For example, pancake-shaped traps 
can produce a two-dimensional (2D) gas, while a cigar-shaped trap can confine it to one 
dimension (ID) [2j. While the static properties of atomic Bose-Einstein condensates 
(BEC) are generally well understood [3] the dynamical behaviour remains an active area 
of study [4]. In the early days of atomic BEC, oscillations of the shape were studied, 
primarily to establish the superfluidity of the condensate, and it was noticed that these 
oscillations were damped I3E] at a rate that depended strongly on the temperature 
fz]- An explanation for this was offered by Landau damping |8l E], in which a low- 
energy excitation of the condensate is dissipated into the thermal cloud by scattering 
phonons from lower to higher energy. Fedichev et al. [mini extended this theory to 
the case of a trapped gas and showed that the damping is determined predominantly 
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by the condensate boundary region, resulting in a different damping rate from that of 
a spatially homogeneous gas. This theory found reasonable agreement with [7], and 
similar agreement was found with the measured damping rate of the scissors mode of 
oscillation [12]. 

Subsequently, Stamper-Kurn et al. ra excited a cigar-shaped condensate to move 
rigidly along its length, out of phase with its thermal component. They saw that 
this second-sound motion [IT] was damped, and noted that collisions neglected in the 
Landau theory might play a role because the hydrodynamicity - the thermal cloud 
collision rate divided by the oscillation frequency - was not small. The damping of 
this mode was also noted in [15] and was studied extensively by Meppelink et al. [16]. 
They found qualitative agreement with [10] at low values of hydrodynamicity, with a 
strongly growing discrepancy at higher values, demonstrating the breakdown of the 
Landau theory at high density. 

Oscillations of long, thin condensates in the ID regime HZ] have very different 
behaviour, with no damping [2] unless corrugation is added to the trapping potential 
[T9] , This raises the question of how the damping evolves from the 3D rate, through 
the crossover regime where no analytic theory currently exists, to a complete absence 
of damping in ID. Oscillation frequencies have been measured in this crossover regime 
[20 , 121] . but not the damping rate. In this article, we measure the damping rate for 
dipole oscillations of a condensate in the crossover regime as a function of temperature, 
and compare our results with measurements of [T6j and the theory of [10J, 11] . We find 
that the oscillations in our experiment persist for longer than expected for a 3D gas and 
propose that this is the consequence of suppressed radial excitations due to the tight 
transverse confinement of the atoms. 

2. Condensate oscillations in a thermal background 

We produce highly elongated, finite temperature condensates with the apparatus 
illustrated in figureH] A magneto-optical trap (MOT) cools and collects 87 Rb atoms a 
few millimetres away from the surface of an atom chip [23]. The MOT is then turned off, 
and the atoms are transferred to a Ioffe-Pritchard trap approximately 110 pm from the 
surface of the chip [2T], 25]. The magnetic trapping field is produced by current in a Z- 
shaped wire on the chip, with its central section along z, together with an external bias 
field along x. The high magnetic field gradient near the centre of the Z-wire gives tight 
radial ( x , y) confinement with a harmonic oscillation frequency of u}p/2n = 1.4 kHz. 
Axial (z) confinement is produced by the currents in the ends of the Z-wire and in the 
end wires (figured]), giving an axial frequency of 3 Hz. 

We cool the trapped gas further by forced evaporation, using an rf field to flip 
the spins of the most energetic atoms so that they are ejected from the trap [26] . By 
sweeping the escape energy down to a few kilohertz above the bottom of the trap, we 
produce an almost pure BEC of approximately 10 4 atoms at a temperature of ~ 150 nK. 
Minor defects in the chip wire cause the current to meander slightly from side to side, 
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Figure 1. Diagram of the apparatus. Cold s 'Rb atoms are delivered from an LVIS 
source [29] to a reflection-magneto-optical trap formed on an atom chip. The atoms 
are passed to a long, thin magnetic trap formed by a current in the Z-wire together 
with a uniform bias field along x. After evaporative cooling, these form a BEC (dark 
blue). A brief ac current in the end wire excites the condensate to oscillate along z, 
as indicated by the green arrows. After some time, the cloud is released and allowed 
to fall for 2 ms under gravity along y, before being imaged along x by absorption of a 
laser beam, shown in red. 


producing small undulations of the trapping potential that make local minima along 
the z axis up to a microkelvin in depth pH EH]. We adjust the centre of the axial trap 
so that the BEC forms in one of these, which is harmonic over a small region, with a 
characteristic frequency of uj z /2h = 10 Hz. The condensed atoms are confined to that 
region, while the higher-energy atoms in the thermal component of the gas explore a 
larger axial range, and experience an anharmonic potential. 

When the rf held is turned off, the atoms warm up at approximately 50 nK/s, 
presumably due to noise in the apparatus. To counteract this, we leave the rf held on, 
so that atoms above some hxed energy are able to leave the trap. Over a few milliseconds 
the cloud comes to equilibrium at the temperature where the heating is balanced by the 
evaporative cooling. We select a desired temperature in the range 150 — 310 nK by 
adjusting the rf frequency. The temperature remains hxed over the next 500 ms, while 
the number of trapped atoms decreases, typically by a few percent. 

Our aim is to observe the oscillation of condensed atoms moving through the 
thermal cloud in order to determine the damping rate as the system equilibrates. To 
resonantly excite axial condensate oscillations, we drive an oscillating current in one of 
the end wires at 10 Hz for two periods. After this time, the condensate’s centre of mass 
is left oscillating with an initial amplitude of ~ 12 pm. The thermal atoms are largely 
unaffected because they explore the region outside the local potential minimum and are 
therefore not resonant with this drive. 

We allow the condensate to oscillate through the thermal background for a time t , 
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Figure 2. Density profiles of ultracold atom clouds and oscillations of the condensate, 
(a) Axial column density profiles measured at three temperatures. Red dots: 
155(3) nK. Blue squares: 251(3) nK. Green triangles: 305(3) nK. Solid lines: fits using 
theory described in the text, which takes into account the irregular potential. Dotted 
lines: profiles of the thermal component of the cloud, determined by the same fit to 
theory, (b-d) Condensate oscillations for the same three temperatures. Points show 
the centre of mass of the condensed component after a period of free oscillation. Lines 
show the fits to the damped sinusoid in ([]]). These fits give (b) 7 = 2.0(6) s -1 at 
155 nK. (c) 7 = 3.8(5) s " 1 at 251 nK. (d) 7 = 5.7(1.2) s " 1 at 305 nK. 


before switching off the trap and imaging the cloud to determine the condensate’s centre 
of mass. By increasing t in 12.5 ms steps over a total of 400 ms, we build up a data set 
of the damped oscillation. We repeat this process for clouds at different temperatures 
which we influence by setting the frequency of the rf field as described above. Thus, we 
observe how the system damps as a function of temperature. 

3. Measuring the temperature, condensate centre of mass, and damping 
rate 

We determine the temperature of the gas, and centre of mass of the condensate from an 
absorption image. To image the atom cloud, we release it from the trap (gravity is up in 
figured]), wait for 2 ms, illuminate it with resonant laser light and view the absorption 
along x using a CCD camera. This image is then integrated over y to obtain the one¬ 
dimensional axial number density profile of the cloud, n{z) = J dxdyn(r). The data 
points in figure[2](a) show axial density profiles measured at three different temperatures. 
At the lowest temperature (red dots), the atoms are nearly all in the condensate, with 
very little signal in the broad thermal background, whereas the profile at the highest 
temperature (green triangles) has a clearly visible thermal population on either side of 
the cloud. 

Our analysis of the cloud profile builds on the method of [30] . The trapping 
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potential is well described by U( r) = ^mu> 2 p 2 + V(z), where p is the radial displacement, 
and V(z) is the potential on axis, including the irregularity caused by the meandering 
current. We determine V{z) from the axial density distribution of cold, non-condensed 
clouds as described in [31]. Knowing U(r), we estimate the number density profile of the 
condensate, n c (r), using the Thomas-Fermi approximation. The profile of the thermal 
component is calculated by integrating the Bose-Einstein distribution over the effective 
potential 2gn c (r ) + f/(r), where the first term is the mean-held energy of thermal atoms 
inside the condensate. The cloud is then allowed to evolve freely for 2 ms to account 
for the period of free fall (though we find that this makes no significant difference to 
the axial profile). We fit this theoretical cloud to the measured density profile n(z) 
in order to determine the temperature, the position of the condensate, and the peak 
condensate number density n c (0). We note that the Thomas Fermi approximation is 
not well satisfied in our 3D/1D condensates, but we find from simulations [32] that 
this method still yields accurate temperatures, while the peak condensate density is 
underestimated, typically by 10%. These fits, shown in hgure[2](a) as solid lines, are in 
excellent agreement with the clouds we observe. For the three clouds that are plotted in 
£gure[2](a), we determined the temperatures 155(3), 251(3) and 305(3) nK. The dotted 
lines show the thermal cloud density within the condensed regions. 

In our experiments, the temperature fluctuates by less than 10 nK from one 
realisation to the next - mainly because of fluctuations in the initial number of 
magnetically trapped atoms - and drifts by less than ±20 nK over an hour. The position 
of the BEC is very stable, fluctuating from shot to shot by less than 1 pm, which we 
associate with mechanical instability of the camera and mirror mounts. It does not drift 
significantly over an hour. 

Figures[2](b-d) show plots of the condensate centre of mass oscillations we measure 
at each of the three temperatures used in hgure[2](a). It takes approximately 30 minutes 
to collect the data points for one plot. We have analysed 33 such time sequences, 
covering a range of temperatures from 150 nK up to 310 nK. In each case, the motion is 
well described by the exponentially damped sinusoid, 

z(t) = kle -7t sin (c at + <f>) + C, (1) 

where A, 7 , oj, 0 and C are fit parameters. Parameters A, 0 and C are independent 
of temperature, and co increases only slightly (by 10 %) over this range of temperatures. 
By contrast, 7 depends significantly on temperature, increasing by a factor of three. 

4. Results and discussion 

In hgure[3]we plot the damping rate 7 as a function of the temperature, each point being 
the result of fitting one oscillation curve. The temperature assigned to one point is the 
mean of the ~ 30 temperatures measured in that curve, and this has a standard error 
smaller than the symbols in the plot. A vertical error bar indicates the ler uncertainty 
in 7 for each fit. We note that the dissipated energy of the oscillation has no significant 
influence on the temperature, because it corresponds to a negligible rise of ~ 2 11 K. 
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T (nK) 

Figure 3. Damping rate measured as a function of temperature for the oscillation 
of our highly elongated BEC. Each point is derived by fitting the oscillation of 32 
cloud images to 0. Vertical error bars show the ler uncertainty in 7 . Horizontal 
error bars are smaller than the symbol. Solid line: a least squares fit of 0 to our 
data gives A v = 3.53(15). Shading indicates the standard error from the fit. Dashed 
line: damping rate given by 0 taking A v = 7, as observed with the 3D condensate of 

mm- 


Damping rates measured in 3D BEC oscillation experiments have generally been 
consistent with the Landau damping theory 0 El 0 m ED CE2J US]. For a 3D trapped 
cloud making small oscillations at a frequency oj close to the trap frequency, this theory 
gives the damping rate as (see equation (18) of Ref. [0]) 


r„ 


A v uj 


k B T 

gV n c( 0 ) 


,3/2 


( 2 ) 


Here A u is a numerical coefficient that depends on which collective mode v is excited, a 
is the s-wave scattering length, n c (0) is the peak number density of the condensate and 
g = 47 t h 2 a/m is the usual nonlinear coupling parameter (the /i of [lOj is the same as our 
gn c (0)). At each temperature our measurements give values for the number density and 
oscillation frequency, from which we construct empirical functions n c (0; T) and uj(T). 
Using these functions, we fit (0) to our data with A v as the only free parameter. The 
result is A v = 3.53(15). The solid line in hgure[3] shows this best fit, with the shaded 
region covering the standard deviation. This theory describes our data well, giving a 
reduced y 2 of 0.90. 

The measurements of Meppelink et al. in [16] involve the same mode as our 
experiment, but the comparison of their damping at low hydrodynamicity with (0) yields 
a coefficient A v = 7 [16_; [33]. The essential difference between these two experiments 
is in the dimensionality of the trapped gas. In [16] the chemical potential was at least 
28 times higher than the radial excitation energy, placing their cloud firmly in the 3D 
regime. By contrast, the chemical potential of our cloud is aproximately twice the 
quantum of radial excitation, which places it in the crossover regime between 3D and 
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hydrodynamicity 


Figure 4. Ratio of measured damping rates 7 to the r„ of ©, taking A v = 7, plotted 
versus hydrodynamicity |33] . Red squares: data of [1~6] . showing an increase in this 
ratio with increasing hydrodynamicity. Blue circles: our data, showing a suppressed 
value of A v and no excess damping up to a hydrodynamicity of 5. 


ID. The temperature dependence of our result indicates that the same Landau damping 
idea still applies, even in this crossover regime, but the density of states, which enters 
through the use of Fermi’s golden rule to obtain (J2J), should be modified to account for 
the quantisation of the radial excitations Physically, the thermal excitations in this 
case are more likely to be along 7 , in which case they cannot contribute to the damping, 
and A u is correspondingly reduced. 

The Utrecht experiment [16] measured the damping over a wide range of 
hydrodynamicity. Following in the spirit of [[16], the red squares in figure® plot the 
ratio of their measured damping rates 7 to the T u of (J 2 ]), with A v = 7 [33], plotted 
versus hydrodynamicity. At low hydrodynamicity, the ratio approaches 1 in their data, 
and 0.5 in our data (blue circles), as discussed above. Further, the Utrecht data shows 
an increase in this ratio as the hydrodynamicity increases, indicating that collisional 
processes, not incorporated in the model of () 2 ]h play an important role in damping this 
dipole mode. Figure® shows that such an increase does not occur in our case. We suggest 
that this too is a consequence of the discrete radial excitation spectrum, which although 
broadened at higher collision rates, remains discrete far above a hydrodynamicity of 5 
and therefore suppresses the ability of thermal-thermal collisions to contribute to the 
damping. 

Following [16], we have taken the hydrodynamicity in figure® to be n t h(v re i)a/uj zi 
where n t h = N t hm 3 ^ 2 uj 2 uj z / (47tA;bF ) 3 ^ 2 is the average thermal atom number density 
experienced by thermal atoms in the harmonic trap, according to the Maxwell 
Boltzmann distribution. The quantity (n re i) = 4 [ksT/{'Krn)^^ 2 is the mean relative 
speed between thermal atoms, and a = 8 tt a 2 is the s-wave scattering cross- 
section. In future, it would be better to derive the thermal density from the 
Bose-Einstein distribution in the harmonic trap, which fixes the mean density at 
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0.55C(|)(27rmfc^T) 3 / 2 /^ 3 , f being the Reimann zeta function. This makes no difference 
to our conclusions here, but will be important for any future quantitative study of the 
corrections to Landau damping. 

In all the damping experiments, the energy in the initial coherent motion is very 
large compared with hco. Indeed, the ratio of these is generally greater than the number 
of atoms in the cloud. It is therefore interesting that the analytical theory reproduces the 
measured 3D damping rates, because the theory assumes a Bogoliubov mode of energy 
hu that is weakly excited. The agreement between experiment and theory indicates that 
the damping rate calculated for weak excitations is still applicable when the excitation 
is strong. 

Collective excitations have been simulated numerically using the method of 
Zaremba, Nikuni and Griffin [34], which makes Hartree-Fock and semi-classical 
approximations to derive a mean field equation for the condensate coupled to a 
Boltzmann equation for the thermal cloud. Simulations by Jackson and Zaremba 
[35, ESI 137] . have proved to be in good agreement with the 3D experiments [123, H, 138] 
respectively. However, in the 3D/1D cross-over where hco p ~ /i ~ /c^T, the quantisation 
of the radial excitations is not well approximated by a semi-classical treatment, as we 
have shown here. A fully quantum treatment may be possible using the perturbative 
approach of [9j[39j, but we are not aware of any such treatment in the 3D/1D crossover 
regime. Our results provide a point of reference for such simulations. In future we 
hope to vary the transverse width of our trap in order to elucidate further the damping 
behaviour in this dimensional crossover region. 
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